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ABSTRACT 

It was recently shown that applying a Gaussianizing transform, such as a logarithm, to the nonlinear 
matter density field extends the range of useful applicability of the power spectrum by a factor of 
a few smaller. Such a transform dramatically reduces nonlinearities in both the covariance and the 
shape of the power spectrum. Here, analyzing Coyote Universe real-space dark matter density fields, 
we investigate the consequences of these transforms for cosmological parameter estimation. The power 
spectrum of the log-density provides the tightest cosmological parameter error bars (marginalized or 
not), giving a factor of 2-3 improvement over the conventional power spectrum in all five parameters 
tested. For the tilt, n^, the improvement reaches a factor of 5. Similar constraints are achieved if 
the log-density power spectrum and conventional power spectrum are analyzed together. Rank-order 
Gaussianization seems just as useful as a log transform to constrain Hs, but not other parameters. 
Dividing the overdensity by its dispersion in few-Mpc cells, while it diagonalizes the covariance matrix, 
does not seem to help with parameter constraints. We also provide a code that emulates these power 
spectra over a range of concordance cosmological models. 

Subject headings: cosmology: theory — cosmology: observations — large-scale structure of universe 
— methods: statistical 



1. INTRODUCTION 

The distribution of matter in the Universe on large 
scales is efficiently quantified by the power spectrum of 
its overdensity fluctuations. This is because to a good ap- 
proximation, the density field is a Gaussian random field, 
at early times (as established by observations of the cos- 
mic microwave background, CMB), or on large scales at 
low redshift. The problem of efficient structure quantifi- 
cation is interesting information-theoretically, and also 
practically, for constraining cosmological parameters. 

While the power spectrum of the overdensity S = 
{p — p)/p is the optimal statistic on the largest scales, on 
translinear scales (0.2 < k/[hMpc^^] < 0.8 at redshift 
z = in a concordance ACDM model), the dark-matter 
density field departs substantially from Gaussianity, and 
the power spectrum covariance i natrix develops a signif- 
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Hen measuring a parameter such as the initial power 
spectrum amplitude, modes in the translinear regime 
are highly correlated, giving little additional constraining 
power when analyzed with larger-scale modes. 

This signals a failure of the power spectrum to describe 
the 6 field fully on these scales, and more practically, 
implies a substantial reduction in its power to constrain 
cosmological parameters. The number of Fourier modes 
grows as for a 3D survey, and so it would be a shame if 
these smaller-scale modes could not be used to constrain 
cosmology. 

Methods have been proposed that reduce the covari- 
ance on translinear scales to varying degrees. These 
include pre- whitening (Hamihon 2000|); remov ing large 
halos from a survey ( |Neyrinck &: Szapudi|2007 ); a Gaus- 



sianizing tr ansform; nonlinear wavelet Wiener filtering 
(Zhang et al. 2011); and div iding 6 by its dispersion in 
few-Mpc cells (Neyrinck 2011). In this paper, by "Gaus- 
sianizing transform" wc mean a function applied to pixels 
in a field (e.g. S measured in cells of some resolution) that 
increases the Gaussianity of 1-point probability density 
function (PDF) of this field. Examples of Gaussianizing 
transforms for the cosmological de nsity field include: a 
logarithmic S — » ln(l-f(5) transform ( Neyrinck et al. 2009 
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2011); rank-order Gaussianization d 
giving an exactly Gaussian distribution by m apping the 
1-point PDF onto a Gauss i an of sonie wid th (Weinberg 



1992 Neyrinck et al. 



2009 



Yu et al.|[2011|; anTa Box 



Cox transform ( |Box &: Cox 1964 " Joachimi et al.||2011 ) 



which can be considered a generalization of the logarith 
mic transform with parameters tunable to give vanishing 
skewness and kurtosis. A related statistic to the rank- 



order-Gaussian ized power spectrum is the copula (jScher- 
rer et al.|[20Tol ) 



I'here are reasons to think that A = l n(l -|- S) would 
be mo re appropriate to analyze than 6. 
^9T] 



Coles fc Jones 

pointed out theoretically that a lognormal PDF 
emerges if peculiar velocities are assumed to grow ac- 
cording to hnear theory. Using a S chrodinger-equation 



framework, Szapudi & Kaiser ( 20d3| found that the vari- 
ance in A is much better-described in tree-level pertur- 
bation theory than the variance in S, suggesting that A is 
closer to linear theo ry. In a study of discreteness effects. 



Romeo et al. (2008) also observed in simulations that the 
first tew moments of A have reduced fractional variance 
compared to 6. 
Going beyond the inherent statistics of A into 



parameter-dependence, Carron (2011) found analytically 
that for a lognormal density field d whose moments de- 
pend on a cosmological parameter, the underlying Gaus- 
sian field A is (often much) more informative about the 
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Fig. 1. — Nonlinear transfer functions of 37 Coyote Universe sim- 
ulations spanning a space of five cosmological parameters. The 
denominators in each panel are no-wiggle linear power spectra, 
without considering effects from the NGP pixel window function. 
The pixel window function causes the small-scale downturns in the 
bottom two panels (and attenuates as well). The dashed blue 
curves show the ratio of to Pnowig • The dotted red curves are 
averages from the 37 simulations. 

parameter t han the lognormal field 5. Also, Joachinii| 
et al. ( [201 1') found that applying a log transform to a 
simulated weak-lensing convergence field allows signifi- 
cantly tighter constraints in a ag-fim parameter space. 
However, they found that adding realistic galaxy shape 
noise in the analysis degrades the constraints both in the 
conventional and transformed convergence fields, reduc- 
ing the gains from Gaussianizing. 

In this paper, we explore the cosmology-constraining 
power of applying three of these transformations to the 
real-space dark-matter overdensity field: a logarithmic 
transform, rank-order Gaussianization, and dividing 5 
by its dispersion in cells. Our analysis ignores the obser- 
vationally relevant issues of shot noise and galaxy bias 
(if the transforms are applied to a gala xy survey), and 



redsh ift-space distortions. In Paper II (Neyrinck et al. 
2011 1, we began the analysis of these issues, to the point 



that we are confident that any cosmological-parameter 
tightening we find in this study will translate to improve- 
ment in a realistic situation as well, although probably 
to a smaller degree. 

2. RESULTS 

In Paper I, we showed that Gaussianizing the low- 
redshift Millennium simulation ( Springel et al. 2005 1 
matter-density field seems to restore a li near shape; here 
we test this a bit more generally. rWang et al. (20111 



studied the scale-dependence of the power spectrum of 
the log-density, Pin(n-5), in renormalized perturbation 
theory. In that paper, in the perturbatively predicted 
P\n{i+S) there is a hint of decreased nonlinearity in its 
shape compared to the conventional power spectrum P5, 
but the perturbative approach does not reach deeply into 
nonlinear scales. 
Fig. [1] shows, for each of 37 high-resolution simula- 



and ^0(15) (the power spectrum of the rank-order- 
Gaussianized field G{5), here mapping 5 onto a Gaussian 
of unit variance) , to the no- wi ggle linear power spectrum 



^'nowig (Eisenstein & Hu 1998j), using that simulation's 



tions in the Coyote Unive rse (Heitmann et al. 2010 
2009[ Lawrence et al. 2010) suite, ratios of Fj, r^\a(i+s) 



cosmology. 

Each curve is, essentially, a nonlinear transfer function, 
for a slightly different cosmology. There are important 
subtleties, though. 

First, the power spectra are not divided by the 
power spectrum of the actual initial conditions, but the 
ensemble-average linear power spectrum (and a no- wiggle 
version of it, at that). Thus, cosmic- variance noise is 
present. However, we dampen this noise, by showing, for 
each density field, the average of the power spectra after 
applying 52 (up to th e first harmonic in e ach direction) 
sinusoidal weig htings ( [Hamilton eral][2006, HRS). 

Second, only the simulation power spectra (the nu- 
merator, but not the denominator) are attenuated by 
nearest-grid-point (NGP) pixel window functions. We 
plot each curve to its Nyquist wavenumber, where the 
attenuation is substantial. We did not correct for the 
pixel window function in the measurements because for 
P\n(i+S) and Pg(5)i changing the resolution does more 
than simply introduce a small-scale attenuation, for ex- 
ample changing the large-scale amplitude. P\n{i+s) a-nd 
Pg{5) do lack the nonlinear upward ramp at fc > 0.1 
Mpc~^ present in P^, but starting at fc w 0.3 Mpc~^, 
they turn down compared to Pnowig- Compare this to 
Fig. 3 of Paper I, in which the denominator is the power 
spectrum of the simulation's exact initial conditions, in- 
cluding NGP attenuation (possible because of the much 
higher mass resolution in the Millennium simulation) . In 
Fig. 3 of Paper I, the ratios do not depart substantially 
from unity. Thus we mainly ascribe the downturns in 
Fig. [T] to the resolution-dependent NGP pixel window 
function. The Ps curves would also be higher than plot- 
ted at high k without the NGP attenuation. 

Finally, P\-a{i+5) and Pg(S)^ generally biased on large 
scales compared to Ps, are multiplied by a factor to line 
up with Pg in the smallest-fc bin shown. For Pg{S)j this 
process is equivalent to setting the variance of the Gaus- 
sian onto which S is mapped so that the large-scale ampli- 
tude of Pg{S) is the same as for Pg. In all further analysis 
below, for simplicity, we use a Gaussian with variance 1. 

We rank-order-Gaussianize 6 before the HRS weight- 
ings are applied, rather than Gaussianizing each weighted 
density field separately. Arguably, it would be fairest 
to Gaussianize each seperately, since the 1-point PDF of 
each weighted density field will slightly fluctuate from the 
PDF of the unweighted field. In real observations, one 
would be confronted with this fluctuated PDF, not the 
PDF from a larger sample. The main reason we Gaus- 
sianize before weighting is that practicically, it is non- 
trivial to Gaussianize the weighted density fleld, which 
has essentially a highly nonuniform (sinusoidal) selec- 
tion function. Our procedure likely slightly overestimates 
the covariance from Gaussianizing the weighted density 
fields seperately, since this would equalize the variances 
in each. This would likely have similar covariance-killing 
effects a s we found divi ding Ps by the density variance 
to have fNeyr inck||2011 1. 

Each simulation, analyzed at redshift 2 = 0, has 
a different set of cosmological parameters, each a 
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plausible (given current observational constraints) con- 
cordance cosmological model. The simulations oc- 
cupy an orthogonal-array-Latin-hypercube in the five- 
dimensional parameter space = flmh'^ € [0.12,0.155], 
ujb = ^bh^ e [0.0215,0.0235], € [0.85,1.05], w e 
[—1.3,-0.7], fJs G [0.61,0.9]. The remaining cosmologi- 
cal parameters, e.g. h, are set to match the tight CMB 
constraint on the ratio of the last-scattering-surface dis- 
tance to the sound-horizon scale. The 1024'^-particle 
simulations have box size 1300 Mpc, fixed in Mpc (not 
h^^yipc) to roughly line up baryon- acoustic-oscillation 
(BAO) features in k among different cosmologies. Their 
resolution is sufficient for power-spectrum measurements 
accurate at sub-percent level at scales down to fc = 
1/iMpc-^ 

All the results in this paper use 256"^ grids, a resolution 
at which shot noise is negligible even for Pin(i+5) and 



G{S)- 



We do not push to smaller scales because even 



at this resolution, 25 of the simulations have handfuls of 
cells with zero particles. Among these 25, the median 
number of zero-particle cells is 19, with maximum 341, 
still <C 256'^. To apply the log transform, we set the 
effective number of particles Nee{N = 0) = 1/2 in zero- 
particle cells, i.e. as though there were half a particle 
in the cell. This equalizes the distance in log-density 
between cells with and 1 particle, and 1 and 2 particles. 
We experimented with changing Ncs{N — 0) by factors 
of two up and down (to 1/4 and 1). Unsurprisingly, this 
had little effect; in the worst-case simulation with 341 
zero-particle cells, Pin(i+s) changed by at most (over all 
k) 0.03%, with typical changes ~0.01%. 

The level of fluctuation in the nonlinear transfer func- 
tion looks substantially greater for Ps than for P\n{i+s) 
and Pg(S) at large k, but around 0.1 Mpc^^, there is 
not much difference. Some of the small-scale fluctuation 
could be from slight cosmology-dependent inaccuracies 
in the no-wiggle transfer function. Note tha t the dashed 
curve, showing the ratio of Pun, from CAMB (Lewis et al. 
2000|, to Pnowig, departs slightly from 1 at large fc. 



By eye, the BAO are of similar amplitude in Pin(i+5) 
and Pg{S) as in Ps- This is not surprising; Gaussianiz- 
ing does not undo the bulk motions that erase small-scale 
BAO wiggles in Ps- One difference, though, is that in Ps, 
the smallest-scale wiggles sit atop the start of the nonlin- 
ear ramp, which suggests that their detection may have 
to compete with a shot-noise-like (on translinear scales) 
one-halo term. Variance in this term can be seen as 
the source of the translinear covariance (jNeyrinck et al.j 
120061). Note that in all panels, the BAO are liliely a bit 
damped or smeared from the averaging over sinusoidal 
weightings, but to a lesser degree than the results be- 
low, which use up to the second harmonic in weightings. 
This dampening should affect all panels equally, and the 
scales of the weightings (just beyond the left edge of the 
plot) are ~ 5 times larger than BAO scales, so these ef- 
fects are probably small. However, we leave a thorough 
quantitative analysis of BAO detection in Gaussianized 
power spectra to future work. 



2.1. Covariance Matrices 

We estimate the cumulative Fisher informat ion (Fisher' 
1935[ [Tegmark et"aL]|1997| [Neyrinck fc Szapudi, 2007) in 



parameters a and /3 over a range of bin indices i G 7?. as 



dp 



where Cn is the square submatrix of the power-spectrum 
covariance matrix C with both indices ranging over TZ. 

EE (AlnP.AlnP,) = (AP.AP,) / (P,) (P^). The in- 
verse of Fai3 then gives the parameter covariance matrix. 

In the bin range TZ, fcmin is the lowest k not di- 
rectly modulated by the sinusoidal weightings (beyond 
-v/S times the wavenumber of the second harmonic), i.e. 
fc,„i„ ^ 27r/(1300 Mpc) x 2^3 = 0.017 Mpc~^ We in- 
vestigate constraints on parameters as femax varies, up to 
the Nyquist frequency, ^ 0.6 Mpc~^. The bins of fcmax 
vary by factors of ^2 (approximately, since each bin's 
fcmax is the mean fcmax in. the bin). 

The power-spectrum covariance matr ices are m e asured 
from the Coyote Universe suite, as in Neyrinck (2011). 
We used the HRS sinusoidal-weightings method, gomg 
up to the second harmonic to get 248 different power 
spectra from each simulation. This gave an estimate of 
the covariance in InP from each simulation. We then 
formed an average covariance matrix across the simula- 
tions, to reduce noise. We averaged the covariances of 
In P instead of P for numerical stability across cosmolo- 
gies; e.g. in linear theory, the covariance in InP does not 
depend on the power spectrum. 

Although using up to second-order weightings and then 
averaging together the covariance estimates among sim- 
ulations beats down the noise in the covariance matrix 
substantially, the noise persists at a level that likely 
somewhat biases our results. This is because noise in 
a matrix that is inv erted generally biases the inverse 
(Hartlap et al. 2007). Ideally, we would correct for the 
noise as suggested in that paper, but a necessary ingre- 
dient, the number of independent samples used for the 
covariance estimate, is not a straightforward quantity in 
the HRS weightings method. Each simulation gives 248 
power spectra of overlapping subsamples; especially at 
low wavenumber, these power spectra are not necessarily 
independent. But this abundance of per haps- redundant 
power spectra has the advantage of reducing the noise to 
a level low enough that e.g. it always provided naively 
invertible covariance matrices. 

Fig. [2] illustrates the effect of noise on the signal-to- 
noise (S/N)^, indicative of the (inverse) effect on param- 
eter constraints, as well. Roughly, (S/N)^ gives the num- 
ber of statistically independent modes in the box, oc fc^ax 
for a Gaussian field. It is measured by setting the deriva- 
tive factors to unity in Eq. ([!]). The faint curves show 
(S/N)^ from each simulation; the bold dashed curves 
show (S/N)^ using the averaged covariance, which we use 

for the results below (proportional to 1/ ^/{S/Wp). Note 
that the scatter in the faint curves is not just from ordi- 
nary cosmic variance, but (predominantly) from the scat- 
ter in cosmological parameters among the simulations. 

Fig. [2] shows that much bias in the Fisher matrix is 
eliminated by going from a single simulation to an av- 
erage over 37. To assess the level of residual bias in 
(S/N)^ after the averaging, we also show an estimate of 
the (S/N)^ for Ps using a noise-free approximate covari- 
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Fig. 2. — The biasing effect of noise in the covariance matrix 
on the Fisher signal-to-noise ratio (S/N)^. The faint curves show 
(S/N)^, for both Pg (black) and Pin(i+s) (green), for each of the 
37 simulations. The bold dashed curves show (S/N)^, using the 
average covariance matrix used for the results below. (S/N)^ for 
Pi5 in the absence of noise may be estimated using the noise-free 
(but approximate) model of |Neyrinck| | |20TT| | (dotted black). The 
proximity of the dashed and dotted black curves suggests i 
level of residual bias in (S/N)^ after averaging. 



nggests a low 



ance matrix. We use the model of Neyrinck (20111, in 
which the covariance on translinear scales comes irom 
scale-independent multiplicative fluctuations. Its only 
ingredient is the pixel-density variance, measurable with 
negligible noise. We averaged together the model covari- 
ance matrices in the same way as the full, measured ones. 
Comparing (S/N)^ for to this model suggests that the 
residual bias is small. Importantly, the bias is also likely 
at the same level for all four power spectra investigated, 
so it probably does not affect our conclusions. Still, we 
keep in mind that our parameter constraints in all cases 
arc likely slightly optimistic. 

Fig. 3 shows the joint correlation matrix, Cij / ^yCuCjj , 
for boffi (P5(/c), Pin(i+5)(A:)) and (Psik), PG(S)ik)). We 
need their cross-covariance matrices below when we an- 
alyze pairs of power spectra together. 

Fig. [4] shows the cross-correlation diagonals, i.e. the 
correlations between Ps{k) and Pin{i+S){k), and Ps{k) 
and PG{s){k), in the same k bins. These cross- 
correlations are related to t he nonlinear propagators 
(Crocce & Scoccimarrol'2006') of each power spectrum, 
since the nonhnear propagator of Ps is unity on linear 
scales. The nonlinear propagator quantifies the mem- 
ory of the particular Fourier phases and amplitudes of 
the initial conditions, as a function of k. For Pg{S), this 
function dips down to ~ 0.2 on translinear scales. 

The way we interpret this figure is that all three power 
spectra have simila r nonlinear propag ators, consistent 



with the findings of Wang et al 



(20111. However, they 



each have different mode-coupling, or, roughly speaking 
"one-halo" terms, correlated to each other only at the 
-20% level. 

2.2. Derivative Terms 

Here we describe how we estimate the derivative terms 
relevant for Fisher analysis, i.e. the vector Di{k) = 
d\nP{k)/dai. The Coyote Universe simulations were 
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Fig. 3. — Joint correlation matrices of {Ps{k) and Pin(i+S){l^))y 
and (Ps{k) and Pg(S){I'^))- The lower-left square is the correlation 
matrix of Ps alone. The upper-right square is the the correla- 
tion matrix of Pin(i+i) alone (upper-left corner) and Pg(S) alone 
(lower-right corner). The upper-left square is the cross-correlation 
between bins of and Phi(i4-a); the lower-right square is the same 
for Pg and Pq(^s)- The cross-correlation diagonals appear in Fig. 
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Fig. 4. — Cross-correlations, measured at the same k bin, between 
Ps{k) and Pin(i+a)(fc), and Ps{k) and Pg(5)(A:). These diagonals 
do dominate their respective cross-correlation matrices at low k. 
At high k though, the curves become small, and are increasingly 
difficult to find amid the off-diagonal covariance. 

set up to span the 5-p arameter space in an optimal way. 
Lawrence et al. (2010) employed sophisticated techniques 
such as principal-component analysis and Gaussian- 
process modeling to produce a precision power-spectrum 
emulator, CosmicEmu. They also used many lower- 
resolution simulations to analyze larger-scale modes, that 
we cannot use here because the Gaussianized power spec- 
tra are more sensitive to particle discreteness. 

Our simpler approach is cruder, but acceptable for our 
purposes. In each k bin, we model the fluctuations away 
from the mean power spectrum of all simulations as a 
linear combination of contributions from each parameter 
fluctuation, i.e. 



lnP"(fc) = (InP(fc)) + D{k)-{a. - a), 



(2) 



where a is a vector in the space of the five parameters 
(InWm, Inwf), n^, w, Inal). a is the mean of a. over 
all simulations. Linear algebra yields an estimate of the 
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Fig. 5. — Power spectra from Coyote Universe simulations 1- 
6 (solid), together with the predictions from Eq. Cosmic- 
variance noise in each case has been dampened suDstantially by 
averaging over 248 HRS sinusoidal weightings. Only the fluctua- 
tion, second term on the RHS of Eq. |2| is shown, i.e. the geometric 
mean among the simulations is divided out. 

derivative terms D{k) from a quintet of only five simu- 
lations, but it is unusably noisy, the signal swamped by 
cosmic variance in each simulation. 

We enhance the signal in two ways. First, in each 
simulation, we use not the raw P{k), but the ensemble- 
average of the power spectra of the density field after 
applying the 248 sinusoidal HRS weightings. This par- 
ticularly squashes fluctuations away from the mean at 
small k, while preserving the overall shape. However, 
the window functions of the weightings likely convolve 
neighboring k modes together somewhat. In particular, 
this probably dampens BAO wiggles a bit for all power 
spectra. 

Second, we estimate D{k) by finding the median in 
each k bin among many estimates of D{k), each mea- 
sured from a quintet sampled from the 37 simulations. 
Although the results visually converge with only ~ 10'^ 
quintets, we use ~ 3 x 10^ quintets (limited by mem- 
ory), looping through all quartets of simulations twice, 
each time choosing a random simulation to complete the 
quintet. 

Fig. [5] shows power spectra from simulations 1-6, com- 
pared with their estimates from Eq. ([2|. Each power 
spectrum is divided by its geometric mean among the 37 
simulations. Typically, the accuracy is at the few-percent 
level, with occasional deviations up to 10%. However, 
these larger deviations could be from cosmic variance in 
the particular simulation. We provide an emulator of 
these four power spectra[^ but caution prospective users 
to note the above caveats. 

Fig. Is] shows D{k) estimated in this way for Pii^ (es- 
timateT using CAMB), P^, Pin(i+5), Pg{S)-: and Ps/o- (the 
power spectrum of i5/(Tceih the 5 divided by its dispersion 
in cells). Except for Pun, and particularly in the latter 
three cases, these curves depend on the resolution used 
(a 256'^ grid here). However, the general trends here 
should hold for different resolutions. For P^, generally 
the results match those measured from the CosmicEmu 
emulator well. However, there is a large amount of noise 
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that produces a discrepancy at small k for D\-au:f,{k). 
This is by far the parameter with the smallest range ex- 
plored; thus, it is not surprising that our rather crude 
method, without further low-resolution simulations for 
low-A: modes, is not quite adequate to explore it. For 
just the parameter w^, we use D\^^i^{k) evaluated from 
CosmicEmu for Ps. To all other power spectra, we add 
the correction i:'£°„''"'''^""(fc) - Anc.,(fc). This seems to 
improve the accuracy (or at least decrease the noise) , but 
we caution that our results for are less accurate than 
for other parameters. 

A few aspects of Fig. [6] should be pointed out. The 
parameter with the most straightforward connection to 
the power-spectrum shape is rig, the tilt. The shape of 
[k) for Pin(i+5) and Pg{S) is nearly as straight as that 
for the linear power spectrum. In contrast, for Pg and 
Pg /cri (k) bends much more substantially at the onset 
of nonlinearity, indicating decreased sensitivity on small 
scales. This supports the claim that Gaussianization and 
the log-transform dramatically reduce nonlinearities in 
the power-spectrum shape. 

A parameter whose simplicity is obscured by unfortu- 
nate notation is In cr| , which we investigate instead of as 
because for all erg, Di^^2[k) = 1 in linear theory. Reas- 
suringly, like CosmicEmu, for this derivative we obtain 
about 1 on linear scales for Ps, as does CosmicEmu. At 
this cell resolution for Pin(i-i-5), this derivative term is 
decreased to about 0.7, indicating decreased sensitivity 
of the mean power spectrum at each wavenumber. It is 
smaller than for Pg because the log transform generally 
decreases the large-scale amplitude, by a factor of about 
exp(CT2^u)lPaper I). 

In Fig. pj all Pg{S) and Pg /a derivative terms are tied 
together at zero at small scales. This is because for 
each, the variance is unity in 1300/256 Mpc « 4 Mpc 
cells. The derivative terms are generally smaller in ab- 
solute value for these power spectra, which translates 
into poorer parameter constraints below than for Ps and 
P\n{i+S), with the exception of the parameter rig. Cu- 
riously, at this resolution, Di^^2[k) is of comparable 
absolute magnitude for Ps/a and Pin(i+5) for small k. 
Naively, one might expect all information about the am- 
plitude to be destroyed in Ps/m in which one divides 
the power spectrum by the variance in grid-cell densities 
(here, of 4 ft," ^ Mpc size, but this holds to some degree for 
8h^^ Mpc cells). However, recall that the amplitude 
is the variance in 8h~^ Mpc volumes in the linearly, not 
nonlinearly, evolved density field; the nonzero _Di„o.2(fc) 
at small scales is apparently from the rise in the nonlinear 
power spectrum in Ps and Ps/a- 

2.3. Error ellipses 

Fig.[7]shows error bars over the set of five cosmological 
parameters, for Ps and P\n{i+s)- The effective volume for 
these results is (1.3 Gpc)V2 « 1.1 Gpc^ « .5 {h-^ Gpcf. 
The factor of two is from the sinusoidal weightings used 
for the covariance matrices, which effectively halve the 
volume. Along the diagonal, the curves are unmarginal- 
ized error bars over single parameters, holding all else 
fixed. Off the diagonal, we examine error bars allowing 
sets of two parameters to vary at a time. The upper plots 
show how error ellipses contract as /cmax increases, while 
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Fig. 6. — Derivative terms Di (k) = 9 In P{k) / dai , quantifying the sensitivity of each power spectrum to the five explored cosmological 
parameters. The left panel is computed from CAMB. The second, S panel shows Di(k) estimated as described in the text (thin), and 
from CosmicEmu (bold). For 5, we use our own estimates of Diik) for further calculations (except for Di^^^, for which we use the 
CosmicEmu estimate). The (bold) curves in the rightmost three panels show our estimates of Di(k), that we use for further results. For 
uif,, a correction of D'^°™^^'^"^ (k) — Di^i^^{k), measured for P^, has been applied. The original ajj measurements, without the correction, 
appear as the thin curves. Note that slight differences between curves are unlikely to be statistically significant, given the occasionally 
visible level of noise. 




Fig. 7. — Error-bar widths and error ellipses for Pj (black) and Pin(i+i) (green), in the Coyote-Universe space of five cosmological 
parameters. Along the diagonal, the curves are unmarginalized (half-)error bars over single parameters, holding all else fixed. Below 
the diagonal, pairs of parameters are considered; errors are shown in each parameter, marginalizing over the other. Solid curves 
correspond to the parameter of a panel's row; dashed curves correspond to the column. Above the diagonal, error ellipses contract as 
fcmax increases. The ellipses, in square panels of side half-length 0.5, are centered at the means of their parameters. An ellipse is drawn 
for each bin of fcmax running from 0.02 to 0.6 Mpc~^, log-spaced by a factor of \/2 (four times sparser than the bins constituting the 
curves). The bold ellipses are at fcmax = 0.1. Outside, analyzing only large scales, Ps and P\-n{i+6) both give similar constraints. Inside 
the bold ellipses, nonlinear scales, where Pin{i+i) excels, are included, up to the innermost ellipse with fcmax ~ 0.6 Mpc~^. 
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the lower plots show how marginalized error bars shrink. 

Constraints obtained analyzing Pin(i+5) are substan- 
tially smaller than for Ps, for all parameters, typically 
by a factor of 2 or 3 if the analysis is pushed to the 
smallest scales shown. The difference is particularly 
large for n^, where the error bar is reduced by a fac- 
tor of 5. Another parameter whose behavior is simple 
to understand is lncr|. As discussed above, Di^^2[k) is 
smaller for P\n(i+s) than for Ps, at all k. Looking at 
the diagonal, unmarginalized plots, this is why the er- 
ror bars are degraded in P\n(i+s) when only linear scales 
are included. However, when pushing into translinear 
scales, the penalty from the decreased derivative term is 
quickly overcome because of drastically reduced cosmic 
variance, resulting in tighter constraints from Pin(i+5) at 
sufficiently small scales. 

Fig.jSjshows the same figure for Pgjfj and Pg(S)- Except 
for the case of the tilt n^, the constraints from Pg{5) are 
weaker than from Pin(i+5), and often even weaker than 
from Pg. This could be surprising given that the covari- 
ance matrix of Pg{S) has the smallest non-Gaussian com- 
ponent, and the highest diagonality, of any of the power 
spectra considered here. The performance of Pg/a is also 
disappointing given the high diagonality of its covariance 
matrix; the performance is also degraded for Pg/^ com- 
pared to Ps for some parameters. For Pg(S) and Ps/ai 
this behavior is from small derivative terms D{k). As 
discussed above, this is largely from the unit variance 
enforced in cell densities for these density fields. 

Other analysis procedures are certainly possible. It 
would be convenient to use Ps on linear scales, but on 
translinear scales, to exploit the reduced nonlinearity in 
the shape and covariance of Pg(S) ■ In principle, this could 
be done by setting the variance of the Gaussian in Pg{S} 
so that the large-scale amplitude of Pg{S) matches that of 
Ps, or equivalently multiplying Pg{S) by a factor to line it 
up with Ps on linear scales. We experimented with this, 
estimating this factor by averaging Ps{k) / PG(5){k) over 
a range of fc, but the uncertainty in this factor produced 
strong covariance on small scales, comparable to that 
of Ps- In fact, these experiments partially motivated 
the fo rm of the Ps covariance matrix found in |Neyrinck| 
( [20lT] ). 

l^ig. [9] shows the results from another possibility, in 
which we analyze Ps and Pin(i+5) together, and Ps and 
Pg(S) together. Generally, the constraints from analyz- 
ing two power spectra together are simply the minimum 
at each fc,nax of the results from analyzing each individ- 
ually. This is unsurprising given the high, but not total, 
degeneracy between the power spectra. However, for the 
combination of Pg and Pg(S)j there are significant gains 
over analyzing each individually for some parameter com- 
binations, but the constraints are never better than for 



shape and in the covariance, degrading parameter con- 
straints. 



P 



ln(l+,5)- 



3. DISCUSSION 



We have explored the sensitivity to cosmological pa- 
rameters of the power spectra of various transformations 
of the overdensity field. Ps, the conventional power spec- 
trum, benefits from being exactly the linear power spec- 
trum on linear scales. Another benefit of Ps is in the sim- 
ple effects of smoothing on it. However, on translinear 
scales it suffers strong nonlinearities, both in the mean 



Pi 



ln(l+5) 



the power spectrum of the log-density, has the 



most cosmology-constraining power of any power spec- 
trum considered here. Typically, pushing to the smallest 
scales analyzed here, constraints in marginalized and un- 
marginalized error bars are a factor of 2-3 smaller than 
for Ps . The generality of this result suggests that it would 
hold for other cosmological parameters as well. This im- 
provement over Ps comes from the high diagonality of 
Pin(i+<5)'s covariance matrix, and from the small depar- 
tures from the shape of the linear power spectrum. In 
particular, the tilt in the linear power spectrum is 
dramatically better-preserved in Pia{i+s) than in Ps, as 
shown in its derivative terms in Fig. [m The log trans- 
form reduces marginalized and unmarginalized error bars 
in Ug by about a factor of 5. 

Pg(S)^ the power spectrum of the rank-order- 
Gaussianized density field, has a covariance matrix even 
more diagonal than P\n(i+s), and exceeds P\n(i+S) in cu- 
mulative signal-to-noise. Also, it can be directl y applied 
in the case of significant discreteness noise (Neyrinck 



et al. 2011). This is unlike Pi 



ln(l+i5)> 



alt hough simple 

modihcatiohs of the log transform are possible to handle 
the problem. Unfortunately for cosmological constraints, 
though, Gaussianization as implemented here enforces a 
unit variance in cell densities. This degrades parame- 
ter constraints, in some cases to levels even worse than 
for Ps- A notable exception is Us, for which constraints 
similar to Pin(i+5) are obtained. 



A p romising approach explored by [Joachimi et al.| 
(20111 employs a Box-Cox transformation, which is a 



generalization of the logarithmic transform that can be 
calibrated to give a distribution with vanishing skewness 
and kurtosis. Perhaps this approach can reduce the non- 
Gaussian covariance to a level similar to the Gaussianiza- 
tion transform, while retaining some information about 
the power spectrum (e.g. its amplitude) on linear scales, 
as in the logarithmic transform. 

Ps/a, the power spectrum of the ratio ^/ctcoII, where 
(Tccii is (5's dispersion in few-Mpc cells, is the final power 
spectrum that we investigate. Ps/a has an impressively 
low non-Gaussian covariance, nearly to the level of Pg(5) ■ 
However, in a similar way as for Pg(S)j dividing by the 
dispersion erases much of the sensitivity to cosmological 
parameters, providing error bars similar to Ps- 

4. CONCLUSION 

We find that applying a nonlinear transform to the 
nonlinear density field can significantly enhance the 
cosmology-constraining power of the power spectrum, 
but apparently only if the transform preserves some 
linear-scale amplitude information. The log transform, 
for example, reduces error bars by a factor of 2-3; for 
the tilt, this factor reaches up to 5. The dramatic reduc- 
tion in nonlinearities in both the power spectrum covari- 
ance (as quantified previously by the cumulative signal- 
to-noise ratio) , and in the power spectrum shape, is what 
accomplishes this. 

In Paper II, we showed that issues from galaxy discrete- 
ness, perhaps the most obvious problem for a logarith- 
mic transform, can be overcome. A modified logarithmic 
transform still enhances the cumulative signal-to-noise 
ratio in the presence of discreteness noise. If the galaxy 




Fig. 8. — Same as Fig.[7] for Pg^^ and Pg(S)- 



sampling is sufficiently dense, the tightening in parame- 
ter constraints found in the current paper will hold when 
applied to observations. 

However, more work is required to investigate the 
cosmology sensitivity of power spectra of Gaussianized 
power spectra in the face of redshift-space distortions 
and galaxy bias. In Paper II, we began this study, but 
more work is required. Generally, fingers of God, present 
in redshift space, smear the density field, reducing the 
non-Gaussianity of the 1-point PDF, and thus somewhat 
decreasing the gains produced by a Gaussianizing trans- 
form. 

Another important issue is whether Gaussianizing 
transforms are of use in detecting BAO. BAG scales 
are only barely translinear (transtranslinear?), even at 
z = 0, so we expect the gains from covariance reduc- 
tion alone to be modest. As we show qualitatively in 
Fig. [l] Gaussianizing transforms do not seem to alter 
the BAO wiggles substantially, with high-order wiggles 
erased similarly as in Pg by large-scale bulk flows. How- 
ever, the smallest-scale wiggle or two that are not washed 
out lie in the regime where the shot-noise-like one-halo 
term is significant, i.e. on the upward ramp in the nonlin- 
ear transfer function for Ps- This suggests that detecting 
the smallest-scale wiggles may be easier for P\n{i+5) and 



^G(S)- 

As one might expect, there are some situations in 
which a log transform could be only marginally useful, 
and some in which it helps substantially. We are not 
aware of a case in which the transform degrades con- 
straints, if the analysis is pushed to sufficiently small 
scales. But even if there is such a case, analyzing the con- 
ventional Pg together with Pin{i+s) or would give tighter 
constraints (if only marginally) than Pg alone. Given the 
simplicity of these transforms, it seems to be well-worth 
using them observationally. This is even in cases that 
we have not directly tested, such as the highly nonlin- 
ear scales of the galaxy power spectrum or correlation 
function, sensitive to galaxy-formation details. In this 
case, one could also try looking at the ratio Pg / Pin{i+S) 
(or even the ratio of the corresponding correlation func- 
tions). 

I thank Istvan Szapudi and Alex Szalay for valuable 
discussions, and Katrin Heitniann and Adrian Pope for 
help accessing the Coyote Universe simulations. I also 
thank the anonymous referee for useful and informative 
suggestions. I am grateful for support from the W.M. 
Keck and the Gordon and Betty Moore Foundations, 
through Alex Szalay. 
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